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Abstract 

We study the edge states for a quantum anomalous Hall system (QAHS) coupled with a spin- 
singlet s-wave superconductor through the proximity effect, and clarify the topological nature of 
them. When we consider a superconducting pair potential induced in the QAHS, there appear 
topological phases with nonzero Chern numbers, i.e., N = 1 and J\f = 2, where Andreev bound 
states appear as chiral Majorana edge modes. We calculate the energy spectrum of the edge modes 
and the resulting local density of states. It is found that the degenerate chiral Majorana edge modes 
for Af = 2 are lifted off by applying Zeeman magnetic field along the parallel to the interface or the 
shift of the chemical potential by doping. The degeneracy of the chiral Majorana edge modes and 
its lifting are explained by two different winding numbers defined at the time-reversal invariant 
point of the edge momentum. 
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I. INTRODUCTION 



Andreev bound state (ABS) has been one of the central issue in superconductivity and 
condensed matter physics. It has been estabhshed that ABS is generated at the edge of 
unconventional superconductor where pair potential changes sign on the Fermi surface 
0]. ABS realized in spin-triplet chiral p-wave superconductor has remarkable properties 
since the ABS has a linear dispersion around ky = 0, where ky is a momentum parallel to 
,|7|]. The time reversal symmetry is broken and spontaneous charge current is 



the surface 



induced along the edge Nowadays, this ABS has been recognized as a chiral Majorana 
edge mode |8|. It is an analogous state to the chiral edge mode of Quantum Hall system 
(QHS) j^. In QHS and its analogous superconducting state, the bulk edge correspondence 
has been discussed based on the Chern number which is one of the topological numbers in 
condensed matter physics flo| . 

It has been found recently that gapped time-reversal invariant systems also can support 
gapless edge states. In HgTe/CdTe quantum well, hehcal edge modes is generated due to 
;he strong spin-orbit coupling. This system is so called Quantum spin Hall system (QSHS) 
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12| . There are also analogous systems in the world of superconductivity. In non- 
centrosymmetric (NCS) superconductors [isl, Q, where the spin-orbit coupling is important 
as in the case of QSHS, the presence of helical Majorana edge modes has been predicted js, 
The helical Majorana edge mode is a special ABSs which is expressed by two counter 



moving chiral Majorana fermions [8[. Instead of charge current, spin current is spontaneously 
generated along the edge. Several new features of spin transport stemming from these helical 



Majorana edge modes have been predicted 



15Nl9|. Furthermore, it has been clarified there 



15Nl9| and without dispersion in 



are several types of helical Majorana edge modes with 

NCS superconductors |2Ch23|| . As well as chiral p-wave superconductors, the gapless edge 
modes in NCS superconductors can be characterized by the bulk topological numbers. 
In general, superconducting states with topologically protected edge states are dubbed as 



topological superconductor 



24 



has been done 



25l | and the classification of topological superconductors 



25|. In particular, chiral Majorana edge modes have been a hot issue in the 



context of topological quantum computing 26|. It has been known that chiral Majorana 



edge mode is generated in the chiral p-wave pair potential 



However, chiral p-wave 



pair potential is easily destroyed by the impurity scattering as compared to s-wave one. 
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Furthermore, the transition temperature of chiral p -wave superconductor Sr2Ru04 



IS 



rather low. Thus, chiral Majorana edge mode generated from simple spin-singlet s-wave 



superconductors is highly desired 29|, |30|]. It has been proposed that a chiral Majorana 



edge mode is produced at the interface of ferromagnet/ spin- singlet s-wave superconductor 
junction on the substrate of three-dimensional topological insulator (TI) 29|]. Also, a simpler 
scheme using the Rashba spin-orbit interaction and the Zeeman field has been proposed |30|. 
The essential point is the simultaneous presence of the broken inversion symmetry by the 
strong spin-orbit coupling and the time reversal symmetry breaking by ferromagnet or the 
Zeeman field. 

There is another way to realize chiral Majorana edge modes by using the chiral edge state 
of QHS attached with spin-singlet s-wave pair potential. Since the external magnetic field 
is very strong, it is difficult to induce pair potential for an ordinary QHS. However, as Qi 



et.al. has proposed 



31| . this difficulty is overcome by considering a quantum anomalous 



Hall system (QAHS), instead of QHS, where the exchange field is not so strong. QAHS 



can be realized by the doping of magnetic impurity in QSHS 32 



In this scheme, the 



presence of the chiral Majorana edge modes can be controlled by the band mass term m, 
chemical potential and the pair potential A. The number of chiral Majorana edge modes 
can be classified by using the Chern number M js]. 

Althogh Qi et. al have proposed an interesting system, the physical properties of Majo- 
rana edge states and their relevance to observables have not been fully clarified yet. For 
this purpose, calculations of local density of states are important since it can be measured 
by angle resolved photoemission spectroscopy or scanning tunneling spectroscopy. In the 
following, we calculate the energy spectrum of the edge states and the resulting local den- 
sity of states for various Chern number A/" in a QAHS coupled with a spin-singlet s-wave 
superconductor. To clarify the difference between the M = 1 and M = 2 states, we apply 
Zeeman magnetic field. We find that when the direction of the magnetic field is parallel to 
the interface, the degeneracy of the two chiral Majorana edge modes in A/" = 2 states are 
lifted off. In order to understand the topological nature in detail, we evaluate the winding 
number of the bulk Hamiltonian. This number corresponds to the number of the zero en- 
ergy state due to the presence of the bulk edge correspondence. We also clarify that the 
degeneracy of chiral Majorana edge modes in A/" = 2 can be lifted off by shift of chemical 
potential from zero. 
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QAH NI ^ 

FIG. 1: (Color online) Schematic illustration of the phase diagram of the QAH+SC hybrid system 



for 1-1=0, proposed by Qi, Hughes and Zhang 



. The X axis denotes the mass term m and the y 



axis denotes the magnitude of the superconducting pair potential A. J\f = 2 and Af = 1 phases 
with A 7^ are topological superconductor (TSC). Af = phase is non topological superconductor 
(NSC). If A=0, the present system changes from QAH state to normal insulator (NI). For nonzero 
A, the energy gap of the bulk closes at two boundaries of the phase diagram. 



The organization of the paper is as follows. In Sec. [TTl we introduce the model of QAHS 
with spin-singlet s-wave superconductor. In Sec. IIII Al , we calculate the energy dispersion 
of the edge state and the corresponding local density of state at = with and without 
Zeeman magnetic field. In Sec. IIII Bl we introduce the winding number at A;^, = in order 
to study the topological property of above edge states. In Sec. IIII C\ we calculate the energy 
dispersion of the edge state for fi ^ 0. In Sec. llVt we summarize the results. 



II. FORMULATION 



We take a simple QAHS in the two-dimensional square lattice model, which is obtained 
jy the replacement p^^y — )■ sinp^ y and — ^ 4 — 2(cospa. + cosp^^) in the model used in 



3l|. Near the F point, these replacements do not change the low energy physical properties 



of the system. Compared to continuum model, square lattice model is convenient when we 
calculate the local density of states. The Hamiltonian describing low energy excitations of 
the quasiparticle is 

y-QAH = ^ij'lhQAH{p)i^p , (1) 
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with 



h-QAH = d{p) ■ (J 

m{p) A(sinpa; — i smpy) 

A{smpx + i sinpy) —m{p) 




(2) 



(3) 



where d{p) = {Asinpx, A sin py,m{p)), cr is Pauh matrix, and m{p) = m + B{A — 2{cospx + 
cos Py)). The band mass term m{p) determines the magnitude of energy shift of up and down 
spin. A, B and m are material parameters. The sign of m changes the topological property 
of the system. Here note that the presence of B term is crucial to create a QAHS. The 
energy dispersion of the above Hamiltonian is symmetric for the mass term m for 5 = 0, 
but is asymmetric for B^{]. In other words, a nonzero value of B makes the sign of m 
meaningful. We take A=B=1 in our following calculation. 

In the following, we consider the proximity effect by an attached spin-singlet s-wave 
superconductor, where pair potential is induced in the QAHS. The system is described by 
the BdG Hamiltonian, 

hQAH{p) - fJ- A 



n 



BdG 



9 Z-^ P 



At 



-P)+ f^, 



(4) 



with 



A 




(5) 



-Pt 

where fi is the chemical potential, A the induced pair potential of the spin-singlet s-wave 
superconductor. 

In order to calculate the local density of states (LDOS) at the edge, we introduce infinite 
potential along the y-ax is as shown in Fig. |2l We calculate the Green function at x = 1 by t- 



matrix formulation 



35| . The system is infinite for the y-direction while it is semi-infinite for 



the x-direction (See Fig. Since translational invariance is absent along the x-direction, 
only the momentum ky along the ^/-direction is a good quantum number. We express the 
Green functions by using spatial coordinate x and x' for fixed ky as follows; 



G{x, x' , ky, u) = Go{x, x' , ky, u) 
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FIG. 2: (Color online) Schematic illustration of QAHS with spin-singlet s-wave superconductor in 
2D system. A red line denotes the infinite potential barrier and a blue line denotes the edge. In 
the actual calculation, the x-direction of the system is restricted finite width with Nx {Nx = 4096). 
We have checked that the finite size effect is negligible. We use periodic boundary condition along 
the y-direction. 



- Go{x,Q,ky,u)^ -Gq{Q,x' ,ky,u) (6) 

Go(0,0, A;y,a;) 

with 

Go{x, x', ky, io) = ^Yl e''^^"""'^Go(A;x, ky, u) , (7) 

Go{h,ky,u) = ^_^'^^^^^^^ (8) 

where is a total number of lattice points for the x-direction. In the right hand side of 
Eq.dnj, the first term denotes the unperturbed bulk Green function, and the second term is 
the scattering effect at the edge. The momentum resolved LDOS at the edge x = 1 N{ky, u) 
is written as 

N{ky,u) = -- (lm[G^^{x,x,ky,u)] + lm[G22{x,x,ky,u)]) (9) 
7r V / x=i 

where the upper suffix R means retarded Green function; replacing u to oj + iS with in- 
finitesimal positive number 6, and the lower indices 11 and 22 indicate the matrix elements. 
And we obtain the LDOS at the edge for energy u measured from the Fermi level as follows 

^H = i^E^(^^'^) (10) 

where Ny is a total number of lattice points for the y-direction. We set A^a;=A^j^=4096 in the 
actual calculations. In the next section, we will show the spectrum N{ky,u) and D{uj). 
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III. RESULTS 



A. edge states for // = 



The minimum value of the bulk energy gap Eg is given at kx = ky = as far as we 



have studied. In this case Eg = \m + A| or \m — A|. We fix A = 0.25 and yU = 0. We 
change the value of m for (a) m = —0.5, (b) m = 0, and (c) m = 0.5 corresponding to 
A/" = 2, A/" = 1 and A/" = phases, respectively. In every case. Eg is 0.25 3l|. The bright 
line between upper and lower bands corresponds to the chiral Majorana edge modes for 
—0.25 < u < 0.25. On the other hand, the background structure with parabolic shape of 
the spectrum with | w |> 0.25 expresses the continuum level originating from energy bands 
of the bulk QAHS. We can see the chiral Majorana edge modes for (a) and (b), but no edge 
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FIG. 3: Momentum resolved LDOS at the edge N{ky,uj) is plotted as a function of uj and ky for 
A = 0.25. (a):m,=— 0.5, (b):m=0, and (c):m=0.5; these three cases are shown on m-A space in 
(d), where (a), (b) and (c) belong to A/" = 2, TV = 1 and M = case, respectively. (d)Schematic 
phase diagram for /i = 0. 



mode in (c) as expected by the Chern number 31 1. 

It is interesting to look at LDOS at the edge D{u) since it can be observed by scanning 
tunneling spectroscopy (STS). In Fig. HI we plot D{uj) with the same parameters used in Fig. 
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FIG. 4: (Color online) LDOS at the edge. (a)m = —0.5, (b)m = 0, and (c)m = 0.5. At these 
parameters, the gap of the bulk bands is 0.25 in all cases. For (a) and (b), there is no energy gap 
in LDOS. 



131 For m = —0.5 with A/" = 2, the resulting D{uj) has a finite value at a; = 0. It has a peak 
at a; = 0. Similar to this case, for m = with A/" = 1, LDOS has a peak structure at a; = 0. 
On the other hand, for m = 0.5, LDOS has a gap structure where D{u) = for | oj |< 0.25. 
The absence of peak structure at a; = is consistent with JV = 0, where there is no chiral 
Majorana edge mode. The presence of the chiral Majorana edge mode seriously changes the 
resulting zero energy LDOS at the edge. Thus, it is possible to check the presence of chiral 
Majorana edge modes by STS. However, there is no qualitative difference between M = 2 
and M = 1 phases as seen from LDOS. From the discussion based on the Chern number, we 
can expect that there are degenerate two edge states in A/" = 2 phase. 

To discriminate A/" = 1 phase from M = 2 phase, we need a further idea: We apply 
the magnetic field to the system to make sure the difference between these two phases. 
We add the Zeeman term ■ a in the original Hamiltonian, where /is is the Bohr 

magneton, H is the Zeeman magnetic field. The magnitude of the bulk energy gap Eg 
is influenced by the Zeeman magnetic field. When H is along y-direction. Eg is given by 
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FIG. 5: (Color online) Momentum resolved LDOS at the edge N{ky,oj) in the presence of magnetic 
field along the y-direction at A=0.25 and fj,B\H\ = 0.15. (a):m=— 0.5 and (b):m=0. 



Eg =1 ^m'^ + {ijb\H\)'^ - A I In FigEl N{ky,uj) is plotted for A = 0.25 and ^ib\H\ = 0.15, 
where magnetic field is applied along the y direction. The resulting Eg is 0.27 and 0.1 for (a) 
m = —0.5 and (b) m = 0, respectively. As seen from Fig. [5](a), for M = 2 phase, a splitting 
of the degenerate two chiral Majorana edge modes appears in the momentum resolved LDOS 
N[ky,uj). On the other hand, for A/" = 1 phase, single Majorana edge mode remains the 
same even in the presence of the magnetic field. Thus, we can discriminate A/" = 1 and 
M = 2 phase through the momentum resolved LDOS at the edge in the presence of Zeeman 
magnetic field. We also found that the splitting of Majorana edge modes for M = 2 phase 
appears only when the direction of the magnetic field is along the y direction. In order to 
understand the orientational dependence of the magnetic field on the energy spectrum, we 
introduce the winding number of the system at A;j^ = in the following subsection. 



B. Winding number 



The winding number W is one of the topological invariants that is well defined for the 



system with odd- number of spatial dimension [36[. It can count the number of zero energy 
state protected by topological property of the Hamiltonian. Although the dimension of our 
present system is two, if we fix ky to a certain value, we can express the present system as 



an effective one-dimensional Hamiltonian 



361 ] ■ In order to define VF, we look for a hermitian 



matrix F which anti-commutes with the Hamiltonian; 



(11) 
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In the presence of the time reversal (TR) symmetry G, and particle-hole symmetry C, we 
can choose F as F = —iCQ. However, in the present system, the time reversal symmetry is 
broken by m{p). Thus, we must find other Fs. Since we found such Fs only for ky = 0,-7r, 
we focus on these cases in the following discussions. The searching of F is done as follows. 
Since F is a 4 x 4 hermitian matrix, it can be expressed by direct products of two Pauli 
matrices, 

Tf,^ = a^®T^ (12) 

where (yU = 0, 1, 2, 3) operates on the spin space, and (/i = 0, 1, 2, 3) on the particle-hole 
space. The suffix indicates the unit matrix, and 1, 2 and 3 the x, y, and z components of 
the Pauli matrices, respectively. In a similar manner, the Hamiltonian can be expressed by 
the basis of F^^. Then from the ( ant i-) commutation relations between the Pauli matrices, 
we find that only F23 and F32 anti-commute with the Hamiltonian. These matrices also anti- 
commute with the Hamiltonian even in the presence of Zeeman magnetic field along the z 
direction since the applied Zeeman term just changes the mass term in the Hamiltonian. 
On the other hand, if we apply the Zeeman magnetic field along the x-axis, only F23 anti- 
commutes with the Hamiltonian. Furthermore, if the direction of the applied magnetic field 
is the y-direction, only F32 anti-commutes. As is shown below, these differences imply that 
the topological nature of the Hamiltonian strongly depends on the direction of the Zeeman 
magnetic field. 

Using these Fs, we can define the winding number for each F. First, we diagonalize these 
F matrices; 

UlTUr=i''^' ' ] (13) 

\ -/2X2/ 

and transform the Hamiltonian to anti-diagonalized form by Ur- 

UlHik)UrJ ''-^'^ (14) 

Then we calculate the determinant of the sub matrix q{k) in Eq.( |T4l) . and denote its real 
part as mi and imaginary part as m2- 

det \q{k)\ = rrii^k) + i'm2{k) (15) 
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Finally, we obtain the winding number W as the following integral; 



W = — 
2tt 



de{k) 



where 



9{k) = argdet \q{k) \ = tan 



^1 m2{k) 



(16) 



(17) 



Here we only consider the winding number at ky 



mi{k) 

since they are found to be zero at 



ky = 71 in the following cases. We denote the winding number for and as W23 and 
W32, respectively. The winding numbers W23 and W32 become visible by plotting mi and 



(a) 



E 



(b) 



E 0- 



-3V=2 
1 




FIG. 6: (Color online) Contour plots of mi{k) and m2{k) for (a) W23 and (b) 1^32. Here we fix ky 
as ky = and changes k^ from — vr to vr. 



7712 in the first Brillouin zone as shown in Fig. El To make it clear, we enlarge the scale of 
Fig. |6] around the origin, in Fig|71 The number of the rotation of the contour around the 
origin just corresponds to the winding number. The obtained 1^23 and 1^32 are summarized 
in Table [H 

From the winding numbers in Table [H we can derive the following results. Let us first 
consider the case without the Zeeman magnetic field. In this case, both of W23 and W32 are 
well-defined. Then, the bulk-edge correspondence tells us that the number of zero energy 
states a.t ky = should be consistent with these winding numbers. As a result, there should 
be two zero energy edge states at /Cy = for A/" = 2 phase, and one zero energy edge state 
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FIG. 7: (Color online) Enlarged plots of Fig. [H 
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TABLE I: Correspondence between the winding number W and the Chern number M. We also 
show the validity of the winding numbers in the presence of the Zeeman magnetic field. W21, can 
not be defined for the applied magnetic field along y-direction, while W32 can not be defined for 
the applied magnetic field along the x-direction. 

for A/" = 1. Here note that = for M = 2 does not necessarily mean no zero energy 
edge states at ky = 0. To be consistent with W23 = 2 for A/" = 2 at the same time, we need 
to have two zero energy edge states aX ky = 0. These two zero modes excellently agree with 
the two- fold degeneracy of the chiral Majorana edge modes found in Figj3]|^a). 
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Now consider the case with a weak Zeeman magnetic field. As was mentioned above, 
W23 becomes ill-defined if we apply the Zeeman magnetic field in the ^/-direction. Thus, the 
number of the zero energy state at ky = is determined solely by VF32 in this case. It is 
found that the remaining winding number W32 takes the same value as that without the 
Zeeman magnetic field, so for A/" = 2 phase, the zero energy states a.t ky = should vanish. 
This result excellently agrees again with the lifting of the degeneracy of the chiral Majorana 
modes in Fig|5t^a). These winding numbers also explains why the degeneracy is not lifted 
off if the Zeeman magnetic field is along the x or z-direction: In these cases, 1^23 remains 
well-defined, so the two degenerate zero modes also remain at ky — 0. 

Finally, we discuss the edge state in A/" = 1 phase. For any weak Zeeman magnetic field, 
at least one of W23 and W32 is well-defined, and both of them take 1. Thus, a single zero 
energy edge mode is always ensured at ky = 0. In this sense, the edge state at ky = in 
J\f=l phase is 'robust' against perturbation. 



C. edge states for finite 

In this subsection, we consider /i 7^ case, which corresponds to a doped QAHS case. It 
has been clarified that the condition for closing of the bulk band gap is as follows 31| : 



+ fx^ = when n^O. (18) 

The momentum resolved LDOS N{ky,uj) is plotted in Figs. M^a) and (b) for (a) /i = 
and (b) fi = 0.3, respectively. (We set here m = 0.3 and A = 0.25.) Since the magnitude 
of the bulk energy gap Eg is given by the minimum value oi Eg =\ m + y^Jj^~+~^ \ and 
Eg =\ m— a//!^ + I, the resulting Eg is 0.05 and 0.09 for (a) and (b), respectively. There 
is no chiral Majorana edge state in Fig. El^a), but it exists in Fig. Mjo). Therefore, with the 
change of ^, the chiral Majorana edge state is generated. In other words, the change of fi 
introduces the quantum phase transition of the topological number from A/" = to A/" = 1. 
Actually, the boundaries of the three phases [i.e. A/" = 0, 1, 2 ) depend on the values of /i in 
addition to m and A. See Fig. [8]^c). 

Now, we focus on the edge state for A/" = 2 with nonzero ^. In Fig. M^a), the momentum 
resolved LDOS Af{ky,u) is plotted for i^l = 0.2, A = 0.25 and m = —0.3. Since the 
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FIG. 8: (Color online) Energy spectrum of the edge modes are plotted as a function of ky with 



m=0.3 and A=0.25 for {a)n = 0, and (b)/i = 0.3. (c)phase 
phases are plotted as a function of A and m for various fi 



aoundary of A/" = 2, TV = 1, and J\f = 



magnitude of the bulk energy gap Eg is expressed as 

Eg is 0.18 in the present case. Comparing Fig. |9]^a) with Fig|3]^a), we find that the degenerate 
two Majorana edge modes reahzed in A/" = 2 phase for /.i = is hfted off by the doping. 

If we change m from m = —0.5 to m = 0, with fixing ^ = 0.2 and A = 0.25, the 
corresponding Chern number changes from A/" = 2 to A/" = 1. In Fig. Mjo), the momentum 
resolved LDOS J^{ky,u) is plotted, where Eg is given by 0.32. The chiral Majorana edge 
mode exists for | a; |< 0.32. By contrast to the edge state of A/" = 2, that of A/" = 1 is robust 
against doping. 



IV. SUMMARY 

In this paper, we have studied the edge states of QAHS coupled with spin-singlet s-wave 
superconductor by the proximity effect. We have calculated the energy spectrum of the 
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FIG. 9: (Color online) Energy spectrum of the edge modes are plotted as a function of ky with 
/i=0.2 and A=0.25 for (a) m = —0.5, and (b) m = 0, respectively, (c) Schematic phase diagram 
for fi ^ 0. The light yellow region corresponds to Af = 2, and the dark red one to Af = 1. The 
dots denotes the parameters corresponding to Figsl9][^a) and (b). 



edge states and the resulting local density of states for various magnitudes of mass term 
m, pair potential A and chemical potential fi. The presence or absence of chiral Majorana 
edge modes influence seriously on the local density of state. We have clarified that it is 
possible to discriminate the state with nonzero Chern number from that with A/" = 
by observing local density of states by STS. Since chiral Majorana edge modes for A/" = 2 
are degenerate each other, it is difficult to discriminate A/" = 1 and M = 2 state by simply 
looking at the energy dispersion relation. To resolve this problem, we apply the Zeeman 
magnetic field. We find when the direction of the magnetic field is parallel to the interface, 
the degeneracy of the present two chiral Majorana edge modes in A/" = 2 states are lifted 
off. Due to the presence of the bulk edge correspondence, the number of the chiral edge 
modes can be evaluated by the winding number defined at ky = 0. For A/" = 1 state, the 
resulting chiral Majorana edge mode are protected by the topological property of the bulk 
Hamiltonian. 
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